Shear viscosity in weakly coupled 7V=4 Super Yang-Mills theory compared to QCD 
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We compute the shear viscosity of weakly coupled J\f=4 supersymmetric Yang-Mills (SYM) the- 
ory. Our result for rj/s, the viscosity to entropy-density ratio, is many times smaller than the 
corresponding weak-coupling result in QCD. This suggests that r\js of QCD near the transition 
point is several times larger than the viscosity bound, rj/s > l/in. 
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The experimental study of heavy ion collisions at the 
Relativistic Heavy Ion Collider has provided remark- 
able and sometimes surprising experimental results. One 
prominent surprise is the large size of radial and elliptic 
flow observed in the collisions 0. Indeed, the strength 
of the elliptic flow observed in non-central heavy ion col- 
lisions is as large as in models where the plasma in the 
early stages after the collision is an ideal fluid 0. (A fluid 
is nearly ideal if local equilibration proceeds on length 
and time scales small compared to the scale of inhomo- 
gcneity of the system.) This is a surprise because the 
plasma droplet created in heavy ion collisions is not much 
larger than its intrinsic microscopic length scale-roughly, 
the de Broglie wavelength of a typical excitation. A 
fluid's nonideality is set by (JintrinsicAnhom.)(?7/s), where 
■q is the shear viscosity and s is the entropy density, (rj/s 
is dimensionless in units with K = 1 = fc B ). There- 
fore, rj/s for the quark-gluon plasma must be numeri- 
cally small to display nearly ideal fluid behavior; a value 
rj/s > 0.2 may be enough to reduce the elliptic flow be- 
low what is observed Q (though the data may actually 
require some non-ideality Q). 

In weakly coupled QCD, which should be valid at suf- 
ficiently high energy densities, rj/s is relatively large, 
rj/s ~ 1/ofg lnctjT 1 - However, the temperatures achieved 
at RHIC (probably below 0.5 GeV) are such that the 
(running) coupling is not small. Attempts to interpolate 
|5j between the behavior of rj/s at high temperature |6| 
and low temperature pj suggest rj/s ~ f at the relevant 
temperatures, which may be too high. Unfortunately, 
the only first-principles technique we have to calculate 
the actual behavior in this strongly coupled region is lat- 
tice QCD, which is fraught with large uncertainties when 
extracting real-time behavior such as rj @- 

Recently, a new theoretical perspective arose. There 
is a theory closely related to QCD, namely 7V=4 super- 
symmetric Yang-Mills theory (SYM), where calculations 
can be performed in the limit of a large number of col- 
ors N c and strong 't Hooft coupling A = g 2 N c by using 
string theory techniques. (The 't Hooft coupling cor- 
rectly accounts for the effective coupling strength given 
the large number of degrees of freedom involved.) The 
shear viscosity of SYM theory for large iV c and A has 
been computed; expressed as the ratio rj/s, it is 



small enough for the fluid to behave almost ideally. 

It has further been conjectured by Kovtun, Son, and 
Starinets, after evaluating rj/s in several related strongly- 
coupled theories, that rj/ s = is in fact a lower bound 
on rj/s in all systems [l(J. Together with the formal simi- 
larities between SYM and QCD, and the expectation that 
QCD is strongly coupled at the temperatures relevant in 
heavy ion collisions, this has led to a belief that QCD 
nearly saturates this bound, and that strongly coupled 
SYM may resemble QCD near the transition and 
may be useful for describing other properties of strongly 
coupled QCD [13, O 13 

This argument is also sup- 
ported by thermodynamic information; lattice calcula- 
tions of the pressure of QCD 0] show that, at a few 
times the transition temperature, the pressure is close to 
| of the value at zero coupling, exactly the ratio obtained 
in strongly coupled SYM theory [T^ . 

To explore whether QCD saturates the viscosity 
bound, we think it is useful to examine more carefully 
how much SYM really behaves like QCD. In particular, 
there is a regime where calculations can be carried out in 
both theories; weak coupling. How close are the values of 
rj/s in SYM theory and in QCD at weak coupling? And 
is the physics which sets the viscosity in the two theories 
the same? This paper will address this question. 

At weak coupling, a gauge theory plasma behaves 
much like a gas of quasiparticles. A rough estimate of 
the viscosity and entropy density are 



"0 ~ Z ra fp v (P+e) ~ l m f p nT , s 



(2) 



with n the number density of excitations and T the tem- 
perature. Hence, rj/s is a measure of the ratio of the mean 
free path for large angle scattering, i m f p , to the thermal 
length 1/T. In practice / m f p is momentum dependent and 
scatterings do not fully randomize a particle's direction. 
To turn the estimate into a calculation, we must use ki- 
netic theory (Boltzmann equations). The single particle 
distribution function /(p, x) evolves according to 



(&+v-son P) x,t) = -c a [f] 



(3) 
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with C a [f] a collision term we explain below. Shear 
viscosity is relevant when the equilibrium distribution 
/o(p,x, t) = (exp[7(_E+u • p)/T] ± 1) _1 varies in space, 
diUj 7^ 0. The stress tensor (in the local fluid frame) 
will be shifted by an amount proportional to this ve- 
locity gradient and 77 is the proportionality constant, 
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Tij—Tij ie(l , = -T](diUj+djUi-^5ijdkUk)-(5ijdkUk- (We 
will not discuss the bulk viscosity £.) The lefthand side 
of Eq. (PJ), at lowest order in diUj, reads 



(d i u j )f (l±f )^ = -C a [f} ) 
while the relevant entry in the stress tensor is 



(4) 



T-—T- 



: ^/(03^f (/-(P,x)-/o( P ,x)). (5) 



Both expressions involve the combination piPj/E. The 
collision term is linear in the departure from equilibrium 
and must be inverted to solve for f — fa- 

Arnold, Moore, and Yaffe have shown how the inver- 
sion of the collision term can be performed in QCD by 
variational methods to determine the viscosity jq. We 
follow their treatment, leaving the full details to the ref- 
erences. For each species, one must model the depar- 
ture from equilibrium by a several parameter Ansatz, 
with trial functions (f^ m '{p). To find the viscosity, one 
needs the integral of p 2 /E against each trial function, 
S am = v a f p {p 2 /E)<j)W(p)f (l±f a ) K the multiplic- 
ity of species a) and the integral moments of the col- 
lision operator, C am . bn = f p v a <f>( m) {p)C a [f b {k) = f + 

fo{l±fo)<f> (n) {k)P 2 (cos6 pk )}, with P 2 (cos0j, fe ) the second 
Legendre polynomial of the angle between p and k. The 
viscosity is p 



>) = 



-l 



15T 



(6) 



where (am) is treated as a single index and C as a matrix. 
Their results in massless QCD, divided by the leading- 
order entropy density s = 2ir 2 T 3 g t /45 (with g* the num- 
ber of bosonic fields plus 7/8 the number of fermionic 
fields, 16+36(7/8) = 47.5 for 3-flavor QCD), are 



.4 



N2g*HB/gy/IQ 



A, B 



34.8. 
46.1. 



4.67 N f = 0, 
4.17 N f = 3. 

(7) 

Massless QCD with Nf — 3 is the case closest to the real 
world, because the temperatures in heavy ion collisions 
have m u , ma, m s < T but m c ^> T. 

Now we apply their technique to Af=A SYM theory. 
This is a theory containing gauge fields for an SU(-/V C ) 
symmetry, 4 adjoint Weyl fermions, and 6 real adjoint 
scalar fields, with Lagrangian density given in [l7| . 

As in QCD, two types of collision processes are rele- 
vant; elastic 2 <-> 2 processes and inelastic effective 1 «-> 2 
processes. The latter are really splitting/joining pro- 
cesses induced by soft scattering in the plasma. 

The vacuum elastic processes have remarkably sim- 
ple squared matrix elements, displayed in Table [I] The 
integral of these matrix elements over possible external 
momenta must be performed by numerical quadratures. 
Those elements containing M t or M u require bosonic, and 
X us and X tu require fermionic, self-energy corrections, 



ss 


^SS 


36(M t + M s + M u ) 


SG 


-» SG 


YIMt 


GG 


-> GG 


4(M t + M s + M u ) 


FS 


-> FS 


48M t + lUXus 


FS 


-> FG 


48X US 


FG 


-> FG 


16M t + 16X US 


FF 


-» FF 


64(M t + Ms + M u ) 



TABLE I: Matrix elements squared, summed on all external 
states of given spin, where S, F, G are spin 0, |, 1 respectively. 

Here M t = 2 [sf+s! + l] , X us = 2 [-£ - J - l] , and M u 

and M s are Mt with u <-> t and s«t respectively. Elements 
for SS -> GG, FF SS, FF SG, and FF -> GG can be 
obtained by crossing. 



but with m 2 -) 



2XT 2 



which are the same as in QCD 
and m| = \T 2 /2. 

The splitting processes are in a way simpler than 
in QCD; all particles are in the adjoint representation 
and all hard particles have the same dispersion relation: 
E 2 = p 2 + XT 2 for p 2 ;» AT 2 . However, there are 
more possible splitting processes, because both gauge and 
Yukawa interactions can induce splitting. The total con- 
tribution to C a m,bn due to splitting processes is 



C a 



,bn 



A 2 



2(2tt)3 



dp 



dk X , 



(8) 
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Z = J ABC {p, k) T(p, k) eV/ T f A ( P )f B (k)f c (p~k) (9) 
x [<Wi m) (p) - SaB^Hk) - 8 a cd> { c\p-k)] 
x [5 bA ^\p) - 5 bB ^\k) - 5 bC ^\p-k)] . 

Here J is the splitting kernel given in Table ITU and T is 
the solution to the integral equation 



T{p,k) = \- 1 
2h - iSEF(h) ■■ 



IW 2h - F < h » 



rf 2 q 



A 2 T 3 



(2tt) 2 q 2 (q 2 +2AT 2 ) 



(10) 
(11) 



5E = 



(F(h+Zq) - F(h) 

l—p,k,p—k 

[p 2 +k 2 +{p-k) 2 ]XT 2 



2pk(p~k) 



Apk(p—k) 



The equation for F accounts for splitting due to multiple 
scattering, with q the transverse momentum exchange 
due to a single scattering and h the non-collinearity be- 
tween p and k, h = p x k. 

Using the procedure of Ref. @ and this collision term, 
we find that the viscosity at next-to-leading log order is 



6.174 



A 2 ln(2.36/VA) ' 



(12) 
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ABC 


V ABC 


3abc (Pj k) 


SFF 


12 


p 2 k(p-k)/p 3 k 3 (p-k) 3 


FSF 


12 


pk 2 {p-k)/p 3 k 3 {p-k) 3 


FFS 


12 


pk(p-kf /p 3 k 3 (p-k) 3 


GSS 


3 


2k 2 (p-k) 2 /p 3 k 3 (p~k) 3 


SGS 


3 


2p 2 {p-k) 2 /p 3 k 3 (p-k) 3 


SSG 


3 


2p 2 k 2 /p 3 k 3 {p~k) 3 


GFF 


4 


k(p-k)(k 2 + (p-k) 2 )/p 3 k 3 (p-k) 3 


FGF 


4 


p(p-k)(p 2 + {p-k) 2 )/p 3 k 3 (p~k) 3 


FFG 


4 


pk(p 2 + k 2 )/p 3 k 3 (p-k) 3 


GGG 


1 


{p 4 +k 4 +(p-k) 4 )/p 3 k 3 (p-k) 3 



TABLE II: Splitting kernels for allowed 3-body processes in 
A/"=4 SYM theory. 
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- Weak coupling 
Weak coupling extrapolated 
strong coupling (leading) 
strong coupling (subleading) 



Coupling A=N c g 2 

FIG. 1: Shear viscosity to entropy density ratio r//s in JV*=4 
supersymmetric Yang-Mills theory (SYM). The dotted curve 
is the weak coupling calculation pushed beyond its likely range 
of validity. 



where as before, A = N c g 2 is the 't Hooft coupling. At 
leading order, rj/s is a complicated function of A which 
must be determined numerically. We resolve the 0(y/X) 
ambiguities in its determination using the procedure of 
Ref. [|J. Our (N c independent!) result is plotted in Fig. 
^ which also shows the strong-coupling asymptotic. The 
dotted part of the weak-coupling curve is where we be- 
lieve that corrections to the weak-coupling calculation 
may exceed the factor-of-2 level, so the curve guides the 
eye rather than being a firm calculation. (In the one 
theory where we have an all-orders calculation of r)/s, 
namely large Nf QCD [lj|, the leading-order and exact 
results deviate by about a factor of 2 when the Debye 
screening mass tod reaches the same value as where we 
switch to a dotted line in Fig. QJ) Similarly, the large- 
coupling asymptotic cannot be trusted where it is not 
close to the large A value of 1/47T. The curves suggest that 
strong coupling behavior sets in around A > 10. Note 
for comparison that the weak-coupling expansion for the 
pressure of SYM theory suggests that it approaches 
the strongly-coupled value at a much smaller value of the 
't Hooft coupling, A ~ 2. 




Coupling A = N c g 2 

FIG. 2: r]/s for SYM theory and for QCD, scaled by the 
dominant A dependence and plotted as a function of A. The 
value in SYM is dramatically smaller than in QCD. 



Our result for weakly coupled SYM theory appears 
rather dramatically smaller than the result for QCD, 
both with and without fermions, at the same coupling, 
as shown very clearly in Fig. [21 Naively, this suggests 
that the viscosity of QCD at strong coupling should be 
of order 7 times larger than that of SYM theory, far from 
the viscosity bound and closer to the values for other flu- 
ids near critical points. However, we should explore this 
conclusion a little more carefully, to try to understand 
how this large difference arose. 

The main physics determining the shear viscosity at 
weak coupling is Coulomb scattering. Neglecting all 
scattering processes but Coulomb scattering changes the 
leading-log coefficient A of Eq. J7J) by less than 3% (0.2%) 
for Nf = 3 QCD (SYM theory). Working beyond loga- 
rithmic order, neglecting all processes but Coulomb scat- 
tering shifts our viscosity result by 0(25%). Therefore, 
to good approximation the physics we must compare be- 
tween theories is the physics of Coulomb scattering. 

Two coupling strengths are relevant in Coulomb scat- 
tering; the coupling of a quasiparticle to gauge bosons, 
and the coupling of that gauge boson to all other de- 
grees of freedom in the plasma. The first coupling 
(summed over available gauge bosons) goes as C^g 2 with 
Cr the relevant group Casimir. In the case of SYM the- 
ory, Cr—Ca—N c ; for QCD it is N c =3 for gluons and 
(N 2 —1)/2N C = | for quarks. The second factor depends 
on the number, representation, and statistics of the other 
degrees of freedom in the plasma, in exactly the combi- 
nation which enters in the Debye screening mass squared. 
Therefore it is natural to expect s/r) ~ C^g 2 {m 2 ^/T 2 ). 

The quarks in SYM theory are adjoint rather than fun- 
damental, leading to about a factor of 2 in the Casimir 
and i in rj/s. But much more importantly, the degree-of- 
freedom count which enters in is substantially larger 
in SYM than in QCD. For instance, for N c = 3, SYM 
theory has 4 x 2 x 8 = 64 fermionic degrees of freedom 
(four Weyl fermion species, consisting of a particle and 
antiparticle in 8 colors) and (6+2) x 8 = 64 bosonic de- 
grees of freedom (6 scalars and 2 gauge boson polariza- 
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0.5 1 1.5 2 2.5 

Screening strength m D /T 

FIG. 3: Leading-order rj/s for SYM and for QCD, scaled 
by CavgAmn/T 2 , which we argue in the text captures the 
dominant coupling dependence, as a function of mu/T. 



tion states times 8 colors), for a total of 128; while in 
Nf = 3 QCD, there are only 3 x 4 x 3 = 36 fermionic 
and 16 bosonic degrees of freedom. The Debye masses 
are correspondingly very different; = 2AT 2 in SYM 
theory, but m 2 , = \\T 2 in A c =3, 2V f =3 QCD. 

Since r\ oc Z m f p which scales as C^ 1 for each species, 
a reasonable way to account for the difference in group 
Casimirs in a theory like iVf=3 QCD, with particles in 
multiple representations, is to take an average of inverse 
Casimirs. Therefore, define 

f~t— 1 _ ^mattorff*matter + C*A ff*adj ,..„.. 
°avg = ~ ' 

9* 

with C ma tter the Casimir for the representation of matter 



fields and <?* m atter, ff*adj the contributions of each type of 
field to defined earlier. Scaling r\ by C^to^/CaT 2 
greatly improves the agreement between different theo- 
ries, as shown in Fig. |3J However, ry/s in SYM remains 
below that in QCD, because SYM theory has interac- 
tions (Yukawa and scalar 4-point) which are absent in 
QCD and which introduce additional scattering chan- 
nels, further lowering r\. While these do not contribute at 
leading-log, they become more important as the coupling 
increases. 

To conclude, weak-coupling comparisons of QCD with 
7V=4 super Yang-Mills theory strongly suggests that 
QCD will not approach the viscosity bound, r\j s~> 1 
close to the QCD phase transition or crossover point. At 
weak coupling, SYM theory has a value of rj/s which is 
about 1/7 that of QCD. This difference arises because 
quarks have a smaller coupling to gluons than the gluon 
self-coupling, since they are in a different group repre- 
sentation; and more significantly, because SYM theory 
has many more degrees of freedom available as scatter- 
ing targets than does QCD. 

Finally, we comment that both the size of rj/s and the 
reliability of thermal perturbation theory seem to rely on 
to d /T 2 rather than on A directly; the combination which 
includes a count of the degrees of freedom in the plasma 
is more relevant to thermal behavior. Relating couplings 
in this way, a s — 0.5 corresponds to A = 4.7, where SYM 
has thermodynamics very close to the strongly-coupled 
behavior but rj/s still nearly 10 times larger. It might 
be reasonable to believe that QCD plasmas relevant at 
RHIC fall in this region. 
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